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Abstract 

A useful semiclassical method to calculate eigenfunctions of the 
Schrodinger equation is the mapping to a well-known ordinary differ- 
ential equation, as for example Airy's equation. In this paper we gen- 
eralize the mapping procedure to the nonlinear Schrodinger equation 
or Gross-Pitaevskii equation describing the macroscopic wave function 
of a Bose-Einstein condensate. The nonlinear Schrodinger equation is 
mapped to the second Painleve equation (Pn), which is one of the 
best-known differential equations with a cubic nonlinearity. A quan- 
tization condition is derived from the connection formulae of these 
functions. Comparison with numerically exact results for a harmonic 
trap demonstrates the benefit of the mapping method. Finally we 
discuss the influence of a shallow periodic potential on bright soliton 
solutions by a mapping to a constant potential. 
PACS: 03.65.Ge, 03.65.Sq, 03.75.-b 

1 Introduction 

In the case of low temperatures, the dynamics of a Bose-Einstein conden- 
sate (BEC) can be described in a mean-field approach by the nonlinear 
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Schrodinger equation (NLSE) or Gross-Pit aevskii equation (see, e.g., [1]) 

1,1 °~ +V(x)+gmx 1 t)\ 2 )^(x,t)=m^-, (1) 



2Mdx 2 v ' / dt 

where g is the nonlinear interaction strength. Stationary nonlinear eigen- 
states of the NLSE satisfying ip(x,t) = exp(—ifit/h)ip(x) fulfill the time- 
independent NLSE 

Analytic solutions of the NLSE are available only for some special cases, 
among these the free NLSE [2,3], arrangements of delta-potentials [4-6] and 
potentials given by Jacobi elliptic functions [7]. Thus there is a great inter- 
est in feasible approximation to the NLSE, among which the most popular 
one is the Thomas-Fermi approximation for the ground state in a trapping 
potential. 

Recently two semiclassical methods for the calculation of nonlinear eigen- 
states were proposed. Mention and discuss other semiclassical approaches to 
the NLSE: Bohr- Sommerf eld quantization [8] and the divergence-free WKB 
method [9, 10] 

In this paper we discuss another semiclassical methods to solve the time- 
independent NLSE approximately based on a mapping to Painleve second 
equation. Special attention will be paid to the shift of the chemical potential 
of bound states due to the nonlinear mean-field energy. This method 
provides good results for excited states and it is fairly easy to understand 
and to use. 

We will focus on the one-dimensional case, which arises, for example, in 
confined geometries (see, e.g., [11] and references therein). Then the effective 
interaction strength is given by g = 2fr 2 aN/Ma 2 L , where a is the s-wave 
scattering length, a± is the transverse extension of the condensate and N is 
the number of atoms [12]. We assume that the wave function is normalized 
as|H| 2 = l. 

For convenience we rescale the NLSE (j2J) such that h = M = l.This yields 
the NLSE in the convenient form 



d 2 ^ 



dx 2 

with q 2 (x) = 2(/i-V(x)). 



q 2 (x)ip(x) + 2g\ip(x)\ 2 ip(x) (3) 
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For the most important case of a harmonic trap V(x) = muj 2 x 2 /2 dis- 
cussed in section 0] this is achieved by the rescaling the variables as 

x' = x/£, ip' = Veip } g' = £Mg/h 2 and//' = /V(M (4) 

with the standard length £ = ^Jh/Muj. The rescaled potential is V(x) = x 2 /2 
and the chemical potential is given in units of hu. To get a feeling for 
the relevant dimensions consider a BEC of 10 4 atoms with transverse width 
a± = lOfim. This yields a scaled nonlinearity of g = +20 for a 87 Rb BEC in 
a trap with axial frequency uj = 2ir x 2Hz and g = —20 for a 7 Li-BEC in a 
trap with uj = 2tt x 100Hz. 

Finally, let us note that the nonlinear Schrodinger equation also describes 
the propagation of electromagnetic waves in nonlinear media (see, e.g., [13], 
ch. 8). 



2 The second Painleve transcendent 

One of the most famous ordinary differential equation with a cubic nonlinear- 
ity is Painleve's second equation or shortly the Pjj equation (see, e.g. [14,15]), 

d 2 q 

^ = 2a0 3 + #, a = ±l. (5) 

The solutions (f>k(%), where the index k refers to the asymptotics at y — > +oo, 
are transcendent. In the linear case, which is found for a = 0, the Pu 
equation reduces to the Airy equation. In fact Airy functions are found in 
the asymptotic limit (see below). 

In textbooks one mostly finds results for the repulsive case a = +1. 
However, asymptotic expansions, which will prove itself as quite useful, are 
available for both cases [16, 17]. For y — > +oo the Painleve transcendent 
4>k{%) vanishes as 

Mv) ~ kM(y) (6) 
We can restrict ourselves to k > 0, since equation (J3J) is invariant under a 
sign change of (j){y). Connection formulae, which relate the asymptotic form 
for y — ► — oo to the form for y — * +oo are well known [16, 17]. For a = +1 
and k > 1 the Painleve transcendent diverges. Otherwise the solution is 
oscillatory for negative y with asymptotics 

MV) ~ %r 1/4 sin Qm 3 / 2 - hd 2 In \y\ - ej + 0{\y\-V A ). (7) 
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Figure 1: The Painleve transcendent 4>k{y) (solid blue lines) for k = 0.5 and 
a = +1 (left) and a = — 1 (right) in comparison with the asymptotic forms 
(JTj) (red dashed line) and kAi(y) (green dash-dotted line). 



The constants depend on the parameter k as 



d\k) 
9{k) 



a 



-ln(l -ak 2 ) 
-ad 2 (k) In 2 + a arg 



r i 



-\d 2 {k) 



TV 

4' 



(8) 
(9) 



The form of the second Painlee transcendent is illustrated in figure ^ We 
plotted the Painleve transcendent 4>k(y) i n comparison with the asymptotic 
expansions (J7|) for y < and © for y > for k = 0.5 and a = ±1. One 
observes that the asymptotic expansions are quite accurate already for small 
values of \y\. 

Furthermore, note that the Pu equation with a = +1 can be written as 
a Hamiltonian system [15] 



d(f) dH dp 
dx dp ' dx 

with the Hamiltonian function 



&H 



(10) 



X 



(11) 
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3 The wedge potential 



As a first illustrative example for the application of the Pjj equation we 
consider the real-valued nonlinear eigenstates in a wedge potential 

V{x) = F\x\. (12) 

This potential might appear a bit artificial, but it provides a natural and eas- 
ily understandable example for the nonlinear quantisation using the Pjj tran- 
scendent. Furthermore, the quantum states of cold neutrons in the earth's 
gravity potential above a hard wall corresponding to a half-wedge were mea- 
sured only recently [18]. 

We consider only real states with a defined parity ip{x) = (—l) n ip{—x), 
such that we can restrict our analysis to the positive real line, x > and 
replace \%jj(x)\ 2 ^(x) by i[>(x) 3 By the means of a scaling y = (2F) 1//3 (x — n/F) 
and ip = |(?| _1 / 2 (2F) 1 / 3 0, the NLSE with the wedge potential is transformed 
to the standard form 

d 2 , 

^ = 2a0 3 + # (13) 

with cr = sign(g). The scaled variable y is negative in the classically allowed 
region Fx < /i such that the wavefunction is oscillatory. In the classically 
forbidden region Fx > one has y > and the wavefunction vanishes as 
4>{y) ~ kAi(y). Note that the differential equation (JT3J) does not depend 
on the nonlinear parameter g explicitly - this dependence is hidden in the 
normalization of <pk{y)- Rescaling the normalization condition H^H 2 = 1 
yields 

+oc \n\ 

\^y)\ 2 ^ = jS^- (i4) 

y(x=0) [*-r ) 1 

The quantisation condition can now be deduced from the asymptotic 
form (J2J) of the Painleve transcendent. Note that the definition of a quantum 
number is not so straightforward as in the linear case, as new nonlinear 
eigenstates can emerge and disappear if the nonlinearity g is changed (see, 
e.g. [19]). However, if we restrict ourselves to the nonlinear eigenstates with 
a linear counterpart and thus a defined parity, the quantum number can be 
identified with the number of zeros of the wavefunction. Thus the relevant 
quantisation condition is that the wavefunction ip n (x) must have n zeros. 
Due to the (anti) symmetry ip(x) = (—l) n ip(—x), the wavefunction assumes 
an extremum (n even) or a zero (n odd) at x — 0. 
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Using the asymptotic form (J7J) of the Painleve transcendent, this condi- 
tion can now be cast into an explicit form. As the asymptotic form of the 
Pjj transcendent is basically given by a sine function, a condition for the 
argument of this sine function directly follows from the conditions on the 
wavefunction. In fact, the argument of the sine at y(x = 0) = —2 1 ^ 3 fi/F 2 ^ 3 
must equal (n + l)vr/2. Inserting this into equation (J7|) yields the relevant 
quantisation condition 

K -j^_ _ -ad (k n ) In [y^j ~ 0{K) = —tt, (15) 

where d{k n ) and 6{k n ) are given by equations (jHJ) and (jHJ), respectively The 
advantage of this method is that the problem of solving a nonlinear boundary 
value problem is reduced to a single algebraic equation. 

However, calculating a nonlinear eigenstate with quantum number n for 
a given value of the nonlinear parameter g is not so easy. In fact one has to 
determine the chemical potential /i so that the quantisation condition (|15j) 
and the normalization condition (|14|) are fulfilled simultaneously. This can be 
achieved by an iterative method. It is much easier, however, to start from a 
fixed value of /i. The quantisation condition ([To]) then yields solutions k n for 
different quantum numbers n. Given these values of k n , one can calculate the 
Painleve functions <pk{y) and the effective nonlinear parameter g n (fi) from the 
normalization integral (|14|) . Rescaling the variables to x and ip again directly 
gives the wavefunction ip{x). 

To test the feasibility of this approach, we consider the nonlinear eigen- 
state n = 4 for a wedge potential with F — 1. The resulting wavefunction is 
shown in figure 121 on the left-hand side (dashed red line) in comparison with 
the numerically exact solution (solid blue line). Both wavefunctions are in- 
distinguishable on the scale of drawing. The right side shows the dependence 
of the chemical potential on the nonlinearity g, again in comparison to the 
numerically exact values. One observes a good agreement. The numerical 
results for the NLSE solutions were obtained using the standard boundary- 
value solver bvp4c of MATLAB. 

The only error in this calculation results from the replacement of the 
Pjj transcendent by its asymptotic form (J7J). This error vanishes rapidly for 
larger quantum numbers n, which is illustrated in figure H3 The extrema of 
the Pjj transcendent are given less accurately by the asymptotic form than 
the zeros. Thus the error is larger for even quantum number n. 
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Figure 2: Nonlinear eigenstate with quantum number n = 6 of the NLSE 
for a wedge potential V(x) = \x\. Left: wavefunction for g = 5, right: De- 
pendence of the chemical potential on the nonlinearities g. The semiclassical 
results (dashed red line) are compared to numerically exact results (solid 
blue line). 

4 The harmonic potential 

Now we want to extend the quantisation method presented in the previous 
section to a more important application - the harmonic trap 



A common method used in semiclassics is a comparison of the Schrodinger 
equation to a well-known differential equation, such as Airy's equation [20, 
21]. Similarly we will map the NLSE for the harmonic trap to the the Pn 
equation. 

We use the mapping ansatz 



well known for the linear Schrodinger equation [21], however with an addi- 
tional scaling constant a. Differentiating twice gives 



We demand that the terms proportional to dcj)/dy cancel, which leads to the 
condition 




(16) 




(17) 
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Figure 3: Solution of the NLSE for a wedge potential: Error of the semi- 
classical calculation of the chemical potential fi sc — \x ex in dependence of the 
quantum number n for g — 5. 



Furthermore the term proportional to d 2 //dx 2 is assumed to be small and 
can be neglected. Substituting the Pu equation © and the NLSE (jHJ) into 
equation (fTSj) finally yields 



-1 



0. 



(20) 



In the linear world, which is given by g = or a = respectively, this 
directly gives a differential equation that determines the mapping y(x). If 
the nonlinear effects are small, we can neglect the nonlinear terms in the 
mapping equation which yields 



dy 
dx 



—q 2 (x) 
y(x) 



(21) 



The scaling constant a is now chosen such that the error due to the neglect 
of the nonlinear terms in the mapping equation (J2~T|) is as small as possible. 
In the fashion of a least squares fit, a is chosen such that the error 



X 



dx / \dx 



4>(y(x)) 6 dx 



(22) 



is minimal. This can be done at the end of the calculation, after <p(y(x)) has 
been determined. 
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Figure 4: Nonlinear eigenstates of the NLSE in a harmonic potential V(x) = 
x 2 /2. Left: Wavefunction of the eigenstate with quantum number n — 10 for 
g — 10, right: Dependence of the chemical potential on the nonlinearities g 
for the eigenstates n = 10 — 16. The semiclassical results (dashed red line) 
are compared to numerically exact results (solid blue line). 



This mapping can now be used to approximately calculate eigenstates 
in symmetric single minimum potentials at x — 0, e.g. a harmonic trap 
V(x) = x 2 /2. For wavefunctions with a linear counterpart, that have a 
defined parity, we can restrict our analysis to x > 0. To avoid a divergence 
at the classical turning point x t the mapping has to be such that y(x t ) = 0. 
Thus the integrating equation ()21)1 yields the mapping in explicit form 



y(x) = ± 



4 / \fWW)W 



2/3 

(23) 



where the — sign is taken in the classically allowed region x < x t and the + 
sign is taken in the classically forbidden region x > Xt- 

The quantisation condition is deduced from the asymptotic form of the 
Painleve transcendent exactly as in section El The only difference is that 
the mapping is now given by equation (J2*3*|) . such that the expression for 
y(x — 0) is a bit more complicated. Thus the relevant quantisation condition 
is given by 

^\y(0)\ 3/2 - lad 2 (k n )\n\y(0)\ - 6{k n ) = (24) 

To test the feasibility of this approach, we consider the nonlinear eigen- 
states for a harmonic potential V(x) = x 2 /2. The result for the eigenfunc- 
tions with n = 10 are shown in figure 0] The left-hand side shows the wave 
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Figure 5: Solution of the NLSE for a harmonic potential: Error of the 
semiclassical calculation of the chemical potential fi sc — fi ex in dependence of 
the quantum number n for g = 1 (blue circles) and g = 10 (red crosses). 



function calculated using the mapping procedure (dashed red line) in com- 
parison with the numerically exact solution (solid blue line). The right-hand 
side shows the dependence of the chemical potential on the nonlinearity g, 
again in comparison to the numerically exact values. One observes a good 
agreement. 

Figure El shows results for different quantum numbers n. The error of the 
semiclassical calculation, i.e. the difference of the semiclassical value for the 
chemical potential \i sc and the numerically exact value fi ex is plotted against 
the quantum number n for g = 1 and g = 10. Except for very small values 
of n and g = 10, for which the reduction to the asymptotic form (JJJ) is not 
valid, one obtains reasonable results for the semiclassical approximation. 

The method introduced above can be extended to asymmetric trapping 
potentials. Then one has to construct solutions around the two classical 
turning points separately, which are matched at a 'mid-phase point' [20]. In 
this spirit the restriction to symmetric or anti-symmetric solutions above is 
nothing but a matching of two solutions at the mid-phase point x = 0. 

In the linear case g = a mapping to a similar potential with two classical 
turning points, in fact the harmonic potential, avoids this matching procedure 
[21]. In the nonlinear case, however, the single-turning point equation Pjj has 
some advantages compared to the NLSE with a harmonic potential because 
the Pi i equation is free of movable branch points and connection formulae 
are well known. 
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Figure 6: Bright soliton solution in a cosine-potential of strength w = —0.2 
for \i = —1. Numerically exact results (solid blue line) are compared to 
results from the mapping technique (dashed red line). 

5 Mapping to a constant potential 

It is well-known that the free NLSE 

( " + 2u(f>{y) - 2g<P(yf = (25) 



dy 2 

has soliton solutions. Bright solitons are found for v < and g < 0, given 

by 

(j>{y) = ^fg sech (V^2T(y - y )) , (26) 

and dark solitons for v < and g > are given by 

<f>(y) = VWg tanh (y/V(y - y Q )) . (27) 

Using a mapping technique as in the previous section we explore the 
effects of a small additional cosine potential on these solitons. In fact we 
consider the NLSE 

+ 2(u - w coa(x)Wx) - 2gip(x) 3 = (28) 

dx 2 

Using again the ansatz (fTTjl and following the lines of reasoning of section |3J 
one arrives at 

[2g<f) 3 - 2v(j>] +2(p-v cos(x))0 - 2ga 2 <p 3 = (29) 
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Figure 7: Shift of the chemical potential of a bright soliton in a cosine- 
potential in dependence of the nonlinearity g e ff for w = —0.2 (left) and in 
dependence of the potential strength w for g e g — 4 (right). Numerically exact 
results (solid blue line) are compared to results from the mapping technique 
(dashed red line). 



Again one chooses the scaling factor a to minimize the difference of the 
nonlinear terms 

,1 2 



X 



ax 



dy 



-i 



(p(y(x)) G dx 



and neglects them in equation (|29|) to arrive at the mapping equation 



dy 
dx 



H — w cos(x) 



V 



(30) 



(31) 



To ensure that the right-hand side is positive, one must always be in the 
classically allowed region (/i > w, thus v > 0) or the classically forbidden 
region (/x < —w, thus v < 0); values in the interval /i G (—w,w) cannot be 
treated within this framework. 

To show the validity of this method we calculate a bright soliton solu- 
tion in a cosine lattice V(x) = wcos(x). Figure |H1 shows the wavefunction 
calculated by the mapping method in comparison to the numerically exact 
solution. One observes a good agreement. 

Furthermore we calculate the dependence of the chemical potential poten- 
tial /j, of such a bright soliton on the nonlinarity g e ff for a fixed value of w and 
on the potential strength w for a fixed nonlinearity g e g. The results obtained 
by the mapping method and the numerically exact results are compared in 
figure One observes a good agreement. 
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